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Abstract 

We study the hadronic production of four-quark states with double heavy 
quarks and double light antiquarks at LHC. The production mechanism is that 
a color anti-triplet diquark cluster consisting of double heavy quarks is formed 
first from the produced double heavy quark-antiquark pairs via gg fusion hard 
process, followed by the fragmentation of the diquark cluster into a four-quark 
(tetraquark) state. Predictions for the production cross sections and their dif- 
ferential distributions are presented. Our results show that it is quite promising 
to discover these tetraquark states in LHC experiments both for large number 
events and for their unique signatures in detectors. 
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1 Introduction 

LHC experiments provide the unique opportunity to explore some exotic states 
of heavy quarks. Among them, particularly interesting states are those four-quark 
states consisting of double heavy quarks and double light antiquarks (or their charge 
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conjugator). The existence of such states can be inferred from the heavy quark 
symmetry. The double heavy quarks in the color anti-triplet state may form a diquark 
cluster by the attractive strong interactions. In the heavy quark limit, the double 
heavy quarks move slowly in a small relative velocity v in the rest frame within a 
smaller distance (1/mv), comparing to the size of light degree of freedom (1/Aqcd)- 
Thus in the tetraquark states the double heavy quarks form a color anti-triplet diquark 
cluster. It contributes the color interactions to double light antiquarks as a color 
source of a heavy antiquark. The other two light antiquarks move around it with 
attractive interactions between them. The dynamics of the light degrees of freedom 
of these tetraquark states are similar with those of the heavy baryons. The picture 
is supported by scrutinizing the typical sizes of real hadrons. The size of the heavy 
quarkonium is around 0.2 — 0.3 fm while that of the light hadrons is around 1 fm. 
The masses of such states can be roughly estimated as the sum of two heavy quarks' 
masses and Aqcd- For the tetraquark states containing the cc, cb, and bb quarks, 
their masses are around 3.4 GeV, 6.8 GeV, and 10.2 GeV, respectively. The flavor 
features of this sort of hadrons are very different from the conventional hadrons. Once 
they are discovered in experiments, it will be an undoubted evidence for the existence 
of the tetraquark states. Some theoretical studies on these states were presented in 

literatureslDElElSlElig. 

The spin of a tetraquark state is a composition of the spins of the four quarks and 
the relative orbital angular momenta between them. For S— waves, all orbital angular 
momenta vanish. Then the spin of the tetraquark state is just the composition of the 
spin of each quark or antiquark. The composition of the spins of double heavy quarks 
may be or 1. However, when two heavy quarks are identical, only spin 1 state 
is allowed due to the antisymmetry by exchanging identical fermions. It is also the 
same case for the light antiquarks sector. In this paper, we are interested in only the 
tetraquark states with all orbital angular momenta vanishing. We denote states by 
TqxQi (^ = 0> 1) 5 where Qi and Q 2 represent their heavy flavor indexes and i is the 
spin of the double heavy quarks subsystem. 

As bound states with quite large masses they are difficult to be produced at usual 
high energy machines. Nevertheless, at LHC, they can be produced efficiently via a 
direct production process which involves effects happened at several distinct distance 
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scales. Firstly, two heavy quark- ant iquark pairs are produced via gluon-gluon fusion 
hard subprocess at distance scale 1/m or shorter. Secondly, for those produced two 
heavy quarks with small relative velocity v there are certain probabilities to form a 
color anti-triplet diquark cluster at distance scale 1/mv. Finally the diquark cluster 
evolves into the tetraquark state via the fragmentation process by picking up two 
light antiquarks from the vacuum at the distance scale 1/Aqcd- 

In this paper, we calculate the hadronic production cross sections of T^ c , T bc (i—0, 
1), and T b \ via gluon-gluon fusion process at LHC Our results show that a number 
of these particles can be produced. We also point out the signatures to detect those 
particles. 

The rest of the paper are organized as follows. In Sec. EJ we present the calculation 
of the subprocess gg — > Tq x q 2 QxQi + X. Sec. [3] devotes to the numerical results of 
the cross sections of the tetraquark states at LHC and at Tevatron. In Sec. SJ we 
give some discussions about the results and the signals in the detectors. 

2 Cross section of gg — > Tq q Q1Q2 + X 

As mentioned above, in the production process of the tetraquark states T Qi q 2 , 
there are three hierarchy distance scales, i.e., 1/m <C 1/mv <C 1/Aqcd- Accordingly, 
the cross sections of the subprocesses for the production of the S— wave tetraquark 
states, -^~Q-i02 ' 

can be factored into three different parts accounting for physical effects 
happening at those distinct distance scales: 

a(gg -+ T QlQa ) = j dU 3 C~ 3 (a s ,Q 1 Q 2 ) |^ 3 (0)| 2 J dxD-^^x), (1) 

where d = 1 for T b ° c and T 6 X C , and d = 2 for T^ b and T^ c ; s is the squared invariant mass 
of the double gluons; dll^ is the Lorentz invariant three-body phase space integral 
element; C^(a s ,QiQ2) is the short- distance coefficient describing the production rate 
of the color anti-triplet point-like Q1Q2 state at the energy scale m or higher; ^3(0) 

is the wave function at the origin of the S— wave diquark state; D% Ti (x) is the 

Q1Q2 

fragmentation function of the diquark into the color-singlet tetraquark state Tq i q 2 . 

The short-distance coefficient can be calculated by perturbative QCD and can 
be expanded in terms of a s at the short- distance energy scale m or higher. The 
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rest parts are non-perturbative effects in nature. To estimate the production cross 
sections, one needs to determine their numerical values. The formation of the diquark 
cluster from the free double heavy quarks is described by the wave function at the 
origin. Their numerical values can be estimated by the potential model. The diquark 
cluster provides an anti-triplet color source as a heavy antiquark. Thus, in the heavy 
quark limit, its fragmentation probability for forming the tetraquark states can then 
be approximately described by that for forming the heavy baryons by a heavy quark. 

In this paper, we calculate the cross sections of hadronic production of the S— wave 
states, T} c) T b \, T b ° c , and T b \, at LHC. We compute the tree level short- distance coeffi- 
cients Cz(a s , Q1Q2) in the leading order a A s using perturbation QCD. By estimating 
the nonperturbative matrix elements, we carry out the numerical calculations of the 
total cross sections. Our results show that a number of these particles can be pro- 
duced. 

We first calculate the short-distance coefficients C^(a s , Q1Q2) at the tree level. 
They are proportional to the squared matrix elements M(gg — > (^1^2)3)- To calcu- 
late the matrix elements, one needs to calculate the subprocess of gg — > Q1Q2Q1Q2, 
with Qi and Q 2 moving in the same 3— velocity and in the color anti-triplet state. At 
the tree level, the production processes involve 36 Feynman diagrams for gg — > bcbc, 
and 72 ones for gg — > bbbb and gg — > cccc. The amplitudes can be classified as 
six gauge-invariant subsets in terms of six independent color bases. Given this, the 
calculations of the amplitudes as well as their squares are straightforward. 

\l / 3(0)'s are the wave functions at the origin of the S-wave diquark clusters. Their 
precise values are difficult to be gained since the large range interaction potential 
between the double heavy quarks in the color anti-triplet state is not very clear while 
the short range one is dominated be Coulomb potential. Reasonably, we take values 
predicted by solving Schodinger equation with Coulomb potential, v (r) = — 2a s /(3r). 
Then the predicted ^3(0) 2 's are 0.143, 0.0382, and 0.0198 GeV 3 for bb, be, and cc 
diquark systems, respectively. Including the confinement part, the wave function 
is squeezed to the central region and hence the wave function at the origin will be 
enhanced. Actually, we have done numerical calculations for the diquark state by 
solving the Schrodinger equation using Coulomb potential plus the linear potential 
fixed in the color-singlet case, the numerical value of wave function at the origin is 
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enhanced by 10% - 30%. More reliable prediction for this nonperturbative number 
can be obtained by some nonperturbative method like lattice QCD. 

We now turn to the fragmentation function of (QiQ2)3-cluster to a tetraquark 
state. As discussed above, the produced heavy diquark cluster in the color anti-triplet 
state provides the same color source as the heavy antiquark to the light antiquarks 
in the limit of the ratio of size of the diquark over that of the light antiquarks in 
the tetraquark state. Thus the fragmentation probabilities to produce the tetraquark 
states Tq i q 2 from the heavy diquarks are the same with that to produce the heavy 
baryons from the heavy quarks. Let's take a QED example to illustrate it. Imagine 
the hydrogen ion or deuterium ion passing material. The probabilities forming the 
hydrogen atom or the deuterium one are the same if the velocities of both ions are 
the same since they possess the same electric charge. The fragmentation function is 
defined in the framework of the infinity momentum where the parton moves in the 
speed of light. 

The fragmentation functions to produce the tetraquark states are nonperturbative 
in nature. Thus the shapes of them can only be described by certain phenomenological 
models [S]. One of the most commonly used models is the Peterson model[7], in 
which the fragmentation function takes the following form: 



N 

Ds -+ r QiQ2 (x) = x[1 _ (1/x) _ W(1 _ x)]2 ( 2 ) 



where £q 1 q 2 is the only parameter determining the shape of the fragmentation func- 
tion; N is the normalization constant. Once the fragmentation probability to produce 
the tetraquark state R is given, N is fixed by the following condition: 

dxD- 3 ^T QiQ2 (x) = R, (3) 

The fragmentation probabilities of c — > A c and b — > A& have been measured in 
e+e" collisions M EE]- According to PDG 2006 [9], R(c -> A c ) = 0.094 ± 0.035, 
and R(b — > A b ) = 0.099 ± 0.017. These results are about 0.1. Therefore as a good 
approximation we may take the fragmentation probability of (QiQ2)3 Tq ± q 2 , the 
value of R in Eq.Q, be 0.1. 

From Eqs.([2]) and ([3]) we know the normalization constant N in the Peterson 
model is dependent on the parameter eg. The model suggested a scaling behavior for 



5 



the parameter €q that is proportional to I/tuq. The e b determined by experiments is 
about 0.003 ~ 0.006 [TUl 12] ■ Using the scaling behavior and taking e b to be 0.004, 
we predict that e bc = (^-fe b ~ 0.0023, e cc = (^-fe b ~ 0.011, e bb = {^?t b = 0.001 
and the corresponding normalization constants are 0.0075, 0.0194, 0.0047. Here we 
take m b = 4.9 GeV, m c = 1.5 GeV and rniQ^) = Tng 1 + rriQ 2 . 

The fragmentation functions are energy scale dependent. Their convolution from 
one energy scale to another satisfies DGLAP equation. However, for the hadronic pro- 
duction of the tetraquark states at LHC, the production is dominated by the smaller 
Pt region of the produced particles, i.e., compatible with the tetraquark masses. For 
smaller Pt production of the tetraquark state, there are no distinguishable energy 
scale differences. Thus the evolution of the fragmentation function used here is not 
significant. Therefore, in the numerical calculations we can take the fragmentation 
functions with the initial energy scale to be the tetraquark mass without evolution. 

3 Numerical results of the cross sections of the 
tetraquark states at LHC 

We now turn to the calculation of the total cross sections with a of the subprocess 
given in Eq. (JTJ). According to the parton model, the cross sections of the processes 
pp — > Tq ± q 2 + Qi + Q 2 + X, can be expressed as: 



/d(T 
dx 1 dx 2 fg 1 (xi 1 fi F )fg a (x2,fMF)^p-{g9 -»■ Tq^Vf) , (4) 



da /",,„, v „ / s da 
dpT~ 

where \xp is the factorization energy scale, Pt is the transverse momentum of the 
Tq i q 2 , and f g (x,fip<) is the distribution function of the gluon in the proton. Here 
we use the cteq61 [13] parton distribution function. There is an uncertainty in the 
calculations arising from the choice of the factorization energy scale //p. Here for 
comparison, we take two different values of //p, i.e., /ip = and /j,r/2, where 
A*r = PtQiQ 2 + m QiQ2 w ith Ptq!Q 2 being the transverse momentum of the diquark. 
With the choices of these factorization energy scales and a s given in [14J for the 
cteq61 [13] parton distribution function, we calculate the total cross sections in CMS, 



ATLAS and LHCb with the invariant mass of the pp system, yfs = 7 TeV and 14 TeV, 
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Table 1: The predicted hadronic production cross sections (in unit nb) of the 
tetraquark states with various Pt cuts at LHC with y/s=14 TeV. The pseudo-rapidity 



cuts 



77 1 < 2.5 for CMS and ATLAS, and 1.9 < rj < 4.9 for LHCb are taken. 







LHC (CMS, ATLAS) 


LHCb 




Tj-CUt 

Pt-cvX 




\r)\ < 2.5 


1.9 < 7] 


< 4.9 


fi F = 




/V2 




/W2 


1 bc 


GeV 


0.460 


0.643 


0.263 


0.374 




5 GeV 


0.167 


0.238 


0.0590 


0.0865 




10 GeV 


0.0286 


0.0429 


0.0077 


0.0118 


1 bc 


GeV 


1.73 


2.42 


1.02 


1.45 




5 GeV 


0.567 


0.822 


0.205 


0.302 




10 GeV 


0.087 


0.131 


0.0235 


0.0356 


1 bb 


GeV 


0.129 


0.193 


0.0697 


0.105 




5 GeV 


0.0758 


0.114 


0.0279 


0.0430 




10 GeV 


0.0220 


0.0339 


0.00598 


0.0093 


T 1 

cc 


GeV 


35.6 


37.8 


24.6 


25.7 




5 GeV 


1.29 


1.84 


0.414 


0.629 




10 GeV 


0.072 


0.113 


0.0178 


0.0274 



respectively. The numerical results of the total production cross sections are listed in 
Tables [1] and [2j The Pt and the rapidity distributions are shown in FigsJUEJ From 
Tables [1] and El we see that the total cross section to produce T} c without Pt— cut 
is much larger than the one with Pt > 5.0 GeV. This can be explained with two 
reasons, i.e., the small peak value in the Pt distribution of the subprocess gg — > T^ c 
and the shape of fragmentation function of cc— cluster. They also result in the lines 
of the P-r— distributions across in Figs.Q, ([2$, (jlj, and (EJ) . 

From Table [TJ we see that for LHC with 300 fb~ l integrated luminosity running 
at 14 TeV, for the production of T^ c , around (3.9 — 5.5) x 10 8 events in CMS and 
ATLAS can be accumulated with kinematic cuts Pt > 5 GeV and \rj\ < 2.5, while this 
number is around (1.2 — 1.9) x 10 8 with kinematic cuts Pt > 5 GeV and 1.9 < 77 < 4.9 
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Table 2: The predicted hadronic production cross sections (in unit nb) of the 
tetraquark states with various Pt cuts at LHC with \/i=7 TeV . The pseudo-rapidity 
cuts \r}\ < 2.5 for CMS and ATLAS, and 1.9 < rj < 4.9 for LHCb are taken. 







LHC (CMS, ATLAS) 


LHCb 




7]-CVX 

P T -CUt 




y < 2.5 

1 / 1 


1.9 < r] < 4.9 












1 bc 


GeV 


0.228 


0.352 


0.117 


0.184 




5 GeV 


0.0766 


0.121 


0.0221 


0.0358 




10 GeV 


0.0119 


0.0195 


0.00243 


0.00407 


1 bc 


GeV 


0.860 


1.33 


0.456 


0.717 




5 GeV 


0.265 


0.416 


0.0776 


0.125 




10 GeV 


0.0364 


0.0599 


0.0075 


0.0126 


1 bb 


GeV 


0.0578 


0.0936 


0.0273 


0.0448 




5 GeV 


0.0328 


0.0535 


0.0098 


0.0164 




10 GeV 


0.0088 


0.0145 


0.00180 


0.00304 


T 1 

cc 


GeV 


20.8 


25.9 


13.3 


16.3 




5 GeV 


0.630 


1.01 


0.166 


0.272 




10 GeV 


0.030 


0.051 


0.00580 


0.0099 
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in LHCb. For the production of T bc , both T bc and T bc need to add together since the 
higher mass state will decay into the lower mass state by emitting a photon. After 
doing this, we see that for 300 fb~ l integrated luminosity running at 14 TeV, for 
the production of T bC) around (2.2 — 3.2) x 10 8 events in CMS and ATLAS can be 
accumulated with kinematic cuts Pt > 5 GeV and \rj\ < 2.5, while this number is 
around (0.8-1.2) x 10 8 with kinematic cuts P T > 5 GeV and 1.9 < r] < 4.9 in LHCb. 
For the production of T bb , around (2.2 — 3.4) x 10 7 events in CMS and ATLAS can 
be accumulated with kinematic cuts Pt > 5 GeV and \r]\ < 2.5, while this number is 
around (0.8-1.3) x 10 7 with kinematic cuts P T > 5 GeV and 1.9 < r] < 4.9 in LHCb. 

From Table [2j we see that for 10 fb" 1 integrated luminosity running at 7 TeV, 
for the production of T^ c , around (0.6 — 1.0) x 10 7 events in CMS and ATLAS can 
be accumulated with kinematic cuts Pt > 5 GeV and \rj\ < 2.5, while this number 
is around (1.6 — 2.7) x 10 6 with kinematic cuts Pt > 5 GeV and 1.9 < r] < 4.9 
in LHCb. For the production of T& c , both T bc and T bc need to add together also. 
After doing this, we see that for 10 integrated luminosity running at 7 TeV, for 
the production of T bc , around (3.4 — 5.4) x 10 6 events in CMS and ATLAS can be 
accumulated with kinematic cuts Pt > 5 GeV and \rj\ < 2.5, while this number is 
around (1.0 — 1.6) x 10 6 with kinematic cuts Pt > 5 GeV and 1.9 < rj < 4.9 in LHCb. 
For the production of T bb , around (3.3 — 5.4) x 10 5 events in CMS and ATLAS can 
be accumulated with kinematic cuts Pt > 5 GeV and \rj\ < 2.5, while this number is 
around (1.0 — 1.6) x 10 5 with kinematic cuts Pt > 5 GeV and 1.9 < r] < 4.9 in LHCb. 

We also calculate the total cross sections at Tevatron with a/s = 1.96 TeV . The 
numerical results of the total production cross sections are listed in Table |3j The 
Pt— distributions and the rapidity distributions are shown in Figs J7] and El 

From Table El we see that for 2 fb~ l integrated luminosity running at 1.96 TeV, 
for the production of T c x c , around (0.8 — 1.6) x 10 5 events at Tevatron can be accu- 
mulated with kinematic cuts Pt > 5 GeV and \y\ < 0.6. For the production of T bc , 
again by adding both T b ° c and T b \ events together, we see that for 2 fb^ 1 integrated 
luminosity, for the production of Ti, c , around (0.4 — 0.8) x 10 5 events at Tevatron can 
be accumulated with kinematic cuts Pt > 5 GeV and \y\ < 0.6. For the production 
of T bb , around 0.4 — 0.7 x 10 4 events at Tevatron can be accumulated with kinematic 
cuts Pt > 5 GeV and \y\ < 0.6. 
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Table 3: The predicted hadronic production cross sections (in unit nb) of the 
tetraquark states with various Pt cuts at Tevatron with y / s=1.96 TeV .The rapidity 
cut \y\ < 0.6. 





p nut 

rf CUt 


T l 

L CC 


L bc 


rpQ 


T 1 
L bb 


Hr 


GeV 


3.01 


0.0832 


0.0219 


0.00436 




5 GeV 


0.041 


0.0168 


0.00487 


0.00182 




10 GeV 


0.00149 


0.00178 


0.00058 


0.00039 




GeV 


5.01 


0.157 


0.0410 


0.00837 




5 GeV 


0.081 


0.0312 


0.0090 


0.00352 




10 GeV 


0.0029 


0.00339 


0.0011 


0.00075 



There are some uncertainties in the predicted numerical results. One arises from 
the ambiguity in choosing the factorization energy scale //p. From the Tables, we see 
that two different values of fip lead to around and less than factor 2 difference in the 
cross section. The next leading order result will reduce this uncertainty. However, it 
will be very hard calculation. Another one arises from the wave function at the origin, 
as discussed in Sec. |2j it will increase the cross sections by 20-70 % including the 
linear confinement potential between the heavy quark. Moreover, all excited heavy 
diquark states will decay into the ground diquark states by emitting photons or 7r's. 
Including those contributions, the production cross sections may be enhanced by 2-3 
times. 

Comparing the Pt distributions at LHC and Tevatron, we see that the differential 
cross sections decrease faster at Tevatron than at LHC with the Pt increase. 

In calculating the accumulated events number, we take the integrated luminosity 
as the corresponding collider running at present or planning luminosity for about one 
year. From the above discussions, we see that the predicted event number of the 
tetraquark at the Tevatron is about 3-4 orders of magnitude lower than that at LHC 
with y/s = 14 TeV for one year. 

The events can be discovered at Tevatron is much smaller than the ones can be 
discovered at LHC, especially when the energy and luminosity of LHC reach the 
maximum. 
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P T (GeV) 



Figure 1: The predicted P T — distributions of the tetraquark states, T^, T} c) T fe ° c and 
T^ c , in CMS and ATLAS with v /s=14 TeV, with the pseudo-rapidity cut \r)\ < 2.5 
and /ii?=/iij/2. 



100 -a 




P T (GeV) 



Figure 2: The predicted P?— distributions of the tetraquark states, T^ b , T} c) T fe ° c and 
TjJ., in LHCb at y/s=lA TeV, with the pseudo-rapidity cut 1.9 < rj < 4.9 and 

/iF=/W 2 - 
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Figure 3: The predicted rapidity distributions of the tetraquark states, T b \, T^ c , T b 
and T b \, with fi F =fi R /2, at LHC with y/s=U TeV. 




P T (GeV) 



Figure 4: The predicted P T — distributions of the tetraquark states, T bb , T* c , T bc and 
T bc , in CMS and ATLAS at ^fs=7 TeV, with the pseudo-rapidity cut \r}\ < 2.5 and 
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10 20 30 

P T (GeV) 



Figure 5: The predicted P T — distributions of the tetraquark states, T^, T} c) T fe ° c and 
TjJ., in LHCb at y^s=7 TeV, with the pseudo-rapidity cut 1.9 < rj < 4.9 and fi F =fi R /2. 




Figure 6: The predicted rapidity distributions of the tetraquark states, T^ b , T} c) T fe ° c 
and with fi F =fi R /2, at LHC with ^/s=7 TeV. 
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Figure 7: The predicted P T — distributions of the tetraquark states, T^, T} c) T fe ° c and 
T^, at Tevatron with ^=1.96 TeV, with the rapidity cut \y\ < 0.6 and fj, F =fi R /2. 




Figure 8: The predicted rapidity distributions of the tetraquark states, T^ b , T} c) T fe ° c 
and Tfc with fi F =fi R /2, at Tevatron with ■ N /i'=1.96 TeV. 
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4 Signature and summary 



The decays of the tetraquark states with double heavy quarks possess unique 
signatures in the detectors. To analyze them in detail is very useful for probing them 
in experiments. For the ground tetraquark states, they can not decay by the strong 
and the electro-magnetic interactions. However, they can decay through the weak 
interaction. The inclusive decay can be regards as the cascade weak decays of the 
heavy quarks while the light antiquarks are spectators. When the leptons can easily 
be identified in detectors, we focus on those semileptonic decays of the heavy quarks 
in the tetraquark states. 

Now we analyze the semileptonic decays of the Tq i q 2 states one by one. By weak 
interaction, the c quark can decay to s quark by emitting the e + or and an 
invisible neutrino. In the T* c , double c quarks decay in this way independently, with 
the two antiquarks as spectators. The emitted same plus sign double leptons and the 
P^-distributions can be used to deduce the background. 

In the T b * c , there are two heavy quarks with different flavors. Both the c and the 
b quarks have the semi-leptonic decay modes. As discussed in the last paragraph, 
the c quark can decay to s quark by emitting, e + or fi + , and an invisible neutrino. 
The b quark can decay to the c quark by emitting the e~, or r~ lepton and the 
corresponding invisible anti-neutrino, followed by the c quark decays as above. As a 
result, the two plus sign leptons and one lepton can be detected from those cascade 
decays, giving the special signature for identifying the T b l c . 

The inclusive semi-leptonic decay of T bb is somewhat more complicated. Each b 
quark can decay to the c quark with a lepton, e~, //"", or r~ and the corresponding 
invisible anti-neutrino, followed by the c quark decays as above. These leptons, two 
leptons and two plus sign leptons, provide the unique signature for identifying the T bb 
state. 

The decay vertexes and the b— tagging method will be very useful for identity 
those tetraquark states. Some other useful decays for identifying those tetraquark 
states are those 3-body or 4-body nonleptonic decays. 

Considering the signatures and the efficiencies of the the detectors, we see that 
LHC with high luminosity running at 14 TeV provides best opportunity to discover 
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doubly heavy tetraquark states. In conclusion, our results show that it is quite promis- 
ing to discover these tetraquark states in LHC experiments both for large number 
events and for their unique signatures in detectors. Once they are discovered in ex- 
periments, it will be an undoubted evidence for the existence of the tetraquark states. 
This will start a new stage in hadron physics. 

Acknowledgments: This work is partly supported by the NSFC with the contract 
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